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Preface 


This  report  differs  greatly  from  the  work  which  I 
had  envisioned  at  the  onset*  The  plan  of  accomplishing 
a rather  complete  analysis  of  plasma  lasers  soon  had  to 
be  shelved*  Unfortunately,  questions  of  plasma  generation 
and  cooling,  material  optimization  and  possible  applica- 
tions are  not  answered  in  this  work*  What  is  determined 
is  thisi  Under  certain  limiting  assumptions,  population 
inversions  will  occur  in  a hydrogen  plasma* 

In  investigating  for  this  report,  I compiled  a rather 
comprehensive  bibliography  of  literature  on  the  plasma 
laser*  Appendix  D is  a cross-reference  of  this  bibliography 
which,  hopefully,  will  facilitate  additional  study. 

At  this  time,  I would  like  to  thank  those  members  of 
the  Physics  Faculty  who,  despite  busy  schedules,  always 
found  time  for  helpful  discussion*  Special  thanks  are 
extended  to  Dr*  D*  Shankland  who  helped  me  through  the  woes 
of  a new  computer  program,  and  to  Mr*  Bruce  Fiene,  FTD/lBTEO , 

0 

who  suggested  this  topic* 

To  my  wife  and  two  daughters,  who  had  to  learn  that 
"Grumpy”  was  not  just  a cartoon  character,  "Thank  you*" 

John  J*  Campbell 
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Abstract 

A rapidly  recombining  plasma  preferentially  popu- 
lates the  upper  energy  levels,  thus  producing  population 

inversions*  Laser  devices  in  which  the  active  medium  is 

y %/  y' 

a recombining  plasma  are  termed  plasma  or  'recombination 
lasers • 

A computer  code  was  developed  to  determine  the  level 
populations  of  hydrogen  as  a function  of  time*  The 
instantaneous  cooling  of  the  free  electrons  in  a completely 
ionized  hydrogen  plasma  was  a „umed*  The  resulting  popu- 
lation inversion  produced  estimates  of  gain  which  were 
sufficient  for  laser  action* 

The  absorption  of  Lyman  radiation  in  the  plasma  was 
also  assumed*  Generally,  the  effect  of  this  absorption 
is  to  decrease  the  duration  and  magnitude  of  the  gain* 
However,  when  the  lower  laser  level  and  the  ground  level 
coincide,  the  resulting  gain  can  be  enhanced  by  this  ab- 
sorption* 

Comparison  of  the  computer  results  shows  reasonable 
agreement  with  other  researchers'  computer  predictions* 

A cross-referenced  bibliography  is  provided  to  assist 
in  further  study  of  plasma  lasers* 
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THE  PLASMA  LASER i 
POPULATION  INVERSION 
IN  HYDROGEN  PLASMAS 

I.  Introduction 

In  1963*  Russian  physicists  Gudzenko*  and  Shelepin 
suggested  that  under  certain  conditions  a rapidly  recom- 
bining plasma  should  act  as  a material  with  a negative 
absorption  coefficient  (Ref  40).  Subsequent  theoretical 
studies 1 conducted  principally  within  the  Soviet  Union, 
expanded  the  uses  of  rapidly  recombining  plasmas  as 
possible  laser  sources.  An  investigation  of  open  litera- 
ture suggested  that  little  work  dealing  with  this  aspect 
of  plasma  physics  was  being  performed  outside  Russian 
boundaries. 

The  purpose  of  this  report,  undertaken  at  the  request 
of  Foreign  Technology  Division,  USAF,  was  to  evaluate  the 
feasibility  of  using  a recombining  plasma  as  a laser  medium. 

Time  constraints  which  were  placed  upon  the  completion 
of  this  report  necessitated  a severe  restriction  of  the 
scope  of  this  study.  A search  of  the  available  open  litera- 
ture and  a computer  analysis  of  a recombining  hydrogen 
plasma  were  performed.  The  resulting  population  inversions 
end  gain  are  the  principal  areas  studied  within  this  report* 

♦Occasionally  transliterated  as  Goodzenko  (Refs  31 1 32 » 

33)  • 

/ 
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Questions  of  actual  plasma  production , cooling  and  external 
coupling  will  not  be  addressed,  except  when  inversion  and 
gain  considerations  require  it*  Discussion  of  complex 
plasmas,  plasma  dynamic  lasers,  and  recombination-chemical 
lasers  will  generally  be  limited  to  providing  available 
references  for  further  study. 

Proceeding  from  an  understanding  of  the  above-mentioned 
limitations,  it  is  possible  to  make  an  assessment  of  lasers 
using  recombining  plasmas* 
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II.  General  Background 


k 


Historical  Background 

In  their  1963  work  and  in  subsequent  publications, 
Gudzenko  and  Shelepin  proposed  that  under  certain  condi- 
tions population  inversions  could  be  achieved  in  rapidly 
recombining  plasmas  (Refs  40,  41,  51)*  The  two  physicists 
maintained  that,  once  this  configuration  was  attained,  the 
plasma  then  should  be  capable  of  serving  as  a laser  medium 
for  incident  electromagnetic  radiation. 

This  recombination  principle,  as  it  is  called,  is 
predicated  upon  the  fact  that  free  electrons  which  are 
rapidly  cooled  in  a highly  ionized  dense  plasma  can  be 
preferentially  captured  into  upper  energy  levels.  In  a 
dense  recombining  plasma,  the  probability  of  collisional 
capture  favors  the  population  of  the  upper  energy  levels 
(Ref  50)*  The  resulting  population  inversion  is  short- 
lived, but  may  be  useable  for  laser  action. 

Subsequent  to  the  formulation  of  the  recombination 
principle  came  the  realization  that  the  numerous  popula- 
tion inversions  might  be  used  selectively  to  control  the 
output  wavelength  of  the  laser  in  different  regions  of  the 
electromagnetic  spectrum.  Theoretical  calculations  indi- 
cated that  laser  action  should  be  attainable  in  the  UV  and 
VUV  regions  where  more  conventional  lasers  did  not  work. 
The  possibility  of  attaining  short  wavelengths  through  the 
application  of  the  recombination  principle  has  spurred  on 
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much  of  the  work  in  this  area  outside  the  Soviet  Union 
(Refs  16,  71 )• 

The  available  Russian  literature  does  not  indicate 
that  recombination  laser  action  is  being  preferentially 
developed  in  any  particular  wavelength  range.  Rather, 
Russian  development  appears  to  be  concentrated  upon  ex- 
tending the  recombination  systems  into  complex  plasmas, 
dynamic  lasers  using  recombination  and  chemical-recombina- 
tion laser s. 

Definition 

Laser  devices  which  use  the  properties  of  a recom- 
bining plasma  are  generally  termed  "plasma  lasers"  or 
"recombination  lasers*"  Unfortunately,  in  the  non-Russian 
literature  and,  indeed,  within  the  Russian  works  the  choice 
of  terminology  tends  to  be  inconsistent  and  may  lead  to 
some  confusion.  The  practice  of  describing  some  gaseous 
laser  systems  as  "plasma  laser"  has  created  some  perplex- 
ity. Within  this  report,  the  terms  "recombination  laser" 
and  "plasma  laser"  will  be  used  interchangeably,  but  only 
to  describe  a system  in  which  the  active  medium  is  a 
recombining  plasma* 

Por  the  purpose  of  eliminating  confusion  between  the 
gas  and  plasma  lasers,  it  is  expedient  to  compare  the  two 
systems  in  general.  In  the  gas  laser,  the  needed  popula- 
tion inversion  results  from  filling  the  upper  laser  levels 
almost  directly  from  the  ground  state  through  collisions 
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with  superheated  electrons*  In  the  plasma  laser*  the 
upper  laser  level  is  filled  from  above  with  electrons 
from  the  continuum.  Within  the  gas  laser*  the  electron 
temperature  (Te)  exceeds  the  ion  temperature  (Tt)  corres- 
ponding to  the  percentage  of  ionization  of  the  gas*  In 
the  plasma  laser*  the  electrons  are  supercooled, 

(Ref  59«848).  The  gas  laser  uses  the  leading  edge  of  the 
ionization  pulse  entering  a rarefied  gas*  while  the  plasma 
laser  works  in  a dense  medium  after  the  ionization  source 
has  been  terminated  (Ref  55*1218).  It  is  impossible  to 
continuously  transition  from  the  gas  laser  scheme  to  the 
plasma  laser  without  destroying  the  gain  characteristics 
of  the  medium*  The  time  immediately  following  the  removal 
of  the  ionization  source  is  termed  the  "afterglow  period" 
(Ref  30)  and  will  enter  later  into  discussions  of  experi- 
mental verification  of  the  theory  of  recombination  lasers* 
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III.  Theoretical  Bftgfcfflrgunfl 


At  this  point i a presentation  of  the  theoretical 
basis  which  is  necessary  for  an  understanding  of  the 
existence  of  population  inversions  will  be  made.  This 
will  be  accomplished  by  describing  the  relevent  processes 
which  occur  within  a dense  optically- thin  hydrogen  plasma. 


Population  inversion 

The  population  inversion  between  the  two  levels  is 


a positive  number, 


AN 


*n  ®m 


(1) 


where  Nn  and  Nm  equal  the  total  populations  (cm“^)  of 
levels  with  principal  quantum  numbers  n and  m (n  > m) i 
and  gj^  and  gm  represent  the  degeneracy  factors  of  levels 
n and  m. 

The  total  populations  Nn  and  Nm  change  as  functions 
of  time.  Consequently,  the  population  inversion  which  is 
available  to  assist  in  the  amplification  of  radiation  may 
change  in  time.  The  complete  time-dependent  solution  for 
the  populations  of  all  the  discrete  levels  within  a plasma, 
even  a "simplified"  hydrogen  plasma,  must  often  rely  upon 
approximations.  The  validity  and  applicability  of  these 
approximations  to  obtain  the  solution  of  a particular  con- 
figuration must  merit  careful  scrutiny. 
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Level  Population 

To  better  understand  the  population  changes  of  a 
discrete  level,  consider  a stationary,  homogeneous, 
hydrogen  plasma  which  is  optically-thin  to  its  own 
radiation*  Let  n and  m be  principal  quantum  numbers 
designating  arbitrary  energy  levels  whose  population 
densities  are  Nn  and  Nm  (cm"3)  respectively*  The  density 
of  free  electrons  will  be  represented  by  Ne  (cm"3),  while 
the  continuum  occupied  by  the  free  electrons  will  be 
indicated  by  c*. 

Using  the  notation  of  Table  1,  the  rate  of  change  of 
any  discrete  level  n may  be  written 

^ * “N  [C  £ V(m,n)  ♦ 2 R(n,m)]  N_ 

dt  n m*n  m<n 


♦ B(n,c*)Ne  ♦ £ A(n,m)} 

m<n 


♦ £ [R(m,n)Ne  ♦ A(m,n)]  NB 

m*n 

♦ £ V(m,n)NeNm  ♦ Ne3  B(c*,n) 
m<n 

♦ A(c*,n)N#2  (2) 

where  a and  n take  on  an  infinite  set  of  values  correspond- 
ing to  the  number  of  energy  levels  within  the  plasma* 

Sq  (2)  must  be  supplemented  by  an  equation  which  expresses 
the  conservation  of  electrons* 

Examination  of  Table  I shows  that  two  additional 
mechanisms  which  affeot  the  population  of  a discrete  level 
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TABLE  I 


I 


have  been  neglected i namely*  photo-ionization  from  the 
level  n,  and  upward  discrete  transitions  through  the  ab- 
sorption of  a photon*  In  the  absense  of  a source  of 
radiation  other  than  the  plasma  itself*  the  neglect  of 
these  terms  is  justifiable  if  the  plasma  is  optically-thin 
(Ref  11 i 155)*  If  the  plasma  does  not  satisfy  this  criter- 
ion, then  the  appropriate  terms  (Refs  68i763-764j  72i450) 
must  be  added  to  Eq  (2)* 

Eq  (2)  may  be  written  in  the  following  matrix  formi 


3t°  * SW»  ♦ Dn  (3) 

which  also  represents  an  infinite  set  of  equations*  The 
ease  of  this  notation,  compared  to  Eq  (2),  is  obvious  and 
it  is  found  frequently  within  the  Russian  literature 
(Refs  14,  89)* 

In  this  formulation,  matrix  (K  } is  termed  the  relaxa- 

nmJ 

tion  matrix  (Ref  54i6o6).  The  diagonal  element,  K^,  is 
equal  to  the  absolute  value  of  the  number  of  transitions 
per  unit  of  time  out  of  level  n to  all  other  possible 
states*  The  non-diagonal  element,  K^,  yields  the  number 
of  transitions  from  state  m into  state  n*  In  terms  of 
the  notation  of  Eq  (2), 

- -£[  £ V(n.m)  ♦ £ R(n,m)]  N, 

m>n  m<n 


♦ B(n,c*)Ne  ♦ £ A(n,m)) 


I 


The  increase  in  the  particles  of  level  n from  the 
continuum  is  described  by  the  term  Dn  where 

Dn  “ Ne3  B<c*»n>  * Ne2  (5) 

The  infinite  sets  of  coupled  differential  equations 
which  are  represented  by  Eqs  (2)  or  (3)  are  unyielding  to 
attempts  at  solution  without  further  simplifications*  The 
next  section  of  this  report  will  discuss  the  approximations 
which  are  frequently  used  to  render  these  equations  manage- 
able* 


Historically,  the  systems  of  equations  which  are 
represented  by  Eqs  (2)  and  (3)  have  been  studied  in  de- 
tail to  characterize  the  decay  of  an  ionized  plasma* 

This  section  will  present  the  approximations  which  are 
specifically  used  for  the  purpose  of  solving  for  the 
populations  of  various  levels*  A reading  of  the  available 
background  literature  on  plasma  lasers  will  disclose  that 
there  is  frequent  mention  of  other  approximations*  In 
Appendix  A,  brief  mention  of  these  various  approaches  will 
be  made  and  references  for  further  study  will  be  presented* 
Limiting  the  Number  of  Equations*  The  number  of 
states  over  which  the  summations  of  Eqs  (2)  and  (3)  must 
be  taken  can  be  reduced  by  recognizing  this  facti  Within 
the  plasma,  there  is  an  effective  lowering  of  the  ionization 
potential  and.  consequently,  above  some  level  rig,  the 
states  are  part  of  the  continuum. 
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Physically,  the  principal  quantum  number  is  the 
value  obtained  from  setting  the  maximum  allowed  electron 
radius  equal  to  the  Debye  screening  length  in  the  plasma* 
Margenau  and  Lewis  (Ref  761 595)  and  Griem  (Ref  39* 141) 
use  for  hydrogen 


*2 


r 1 i*/» 

*e 

4nNTe2a£ 


(6) 


where  aQ  equals  the  Bohr  radius,  and  NT  equals  the  total 
particle  density  (cnf^). 

Additional  simplification  may  be  introduced  by  using 
the  work  of  Hinnov  and  Kirschberg  (Ref  64 1 798) • There 
exists  a discrete  level,  n^,  above  which  all  the  discrete 
excited  levels  are  in  Saha  equilibrium  with  the  free 
electrons.  Level  has  this  property 1 The  number  of 
electrons  which  will  recombine  when  transferred  down 
through  it  from  the  continuum  is  equal  to  the  number  of 
electrons  which  will  eventually  ionise  when  passing  up 
through  it  from  lower  levels*  The  principle  of  detail 
balance  implies  that 

-ws  - V‘)A*]  ’ 1 ( 


where  B_ 
*1 


and  E 

nl 


1 


are  the  energies  of  levels  n^  and 


nj-1  respectively* 

Iq  (7)  may  be  solved  by  substituting  for  the 
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* 


degeneracy  factors  of  hydrogen  and  assuming  that  n^x> 1 • 
Expansion  and  algebraic  manipulation  yield 


where  Ry  equals  Rydberg's  constant  (eV),  and  Te  represents 
the  temperature  of  the  electrons  (eV) • 

The  discrete  levels  above  n^  but  below  are  assumed 
to  obey  Saha's  distribution  function 

Nn  * n2  Ne2(h2/2n,neTe^^2  ®*P(lEnlAe)  (9) 

where  h is  Planck's  constant.  me  is  the  mass  of  the 
electron*  and  En  is  the  energy  of  level  n* 

The  net  result,  then,  is  to  solve  Eqs  (2)  or  (3) 
over  a limited  number  of  discrete  states,  nj,  above  which 
there  is  a quasi-continuum  of  discrete  states,  (rig  - n^). 
lying  below  the  actual  continuum*  It  should  be  noted  that 
Eq  (8)  is  not  always  used  to  establish  the  value  of  n^* 
Frequently,  an  arbitrary  selection  is  made,  the  calculations 
performed,  and  the  answers  compared  to  solutions  correspond- 
ing  to  another  choice  of  n^*  The  n^  which  is  finally  chosen 
is  the  smallest  value  that  does  not  significantly  alter 
the  results  (Refs  59»  79) • 

Change  Aq  Electron  Temperature  &£  & Function  Time ■ 

As  previously  noted,  many  of  the  coefficients  which  appear 
in  Eqs  (2)  and  (3)  are  functions  of  the  electron  tempera- 
ture* During  the  rapid  oooling  of  the  plasma,  it  is 
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reasonable  to  expect  that  the  time  history  of  the  electrons 
will  reflect  not  only  the  effect  of  the  cooling  mechanism 
chosen(  but  also  the  heating  effects  of  a recombining 
plasma#  To  date,  attempts  to  allow  for  temperature  change 
in  time  have  fallen  into  two  categories#  The  first  con- 
siders the  temperature  of  the  electrons  as  constant  at 
some  arbitrary  final  value  during  the  entire  recombination 
process t the  other  allows  for  temperature  variation  from 
plasma  processes  by  adding  to  Eqs  (2)  and  (3)  an  equation 
for  the  time  rate  of  change  of  the  electron  temperature 
when  neutral  heavy  particles  are  introduced  as  a thermo- 
static gas  (Ref  94i1062)#  The  former  approach  obviously 
simplifies  the  mathematics  of  the  problem  and  will  be  used 
subsequently  in  this  report#  The  impact  of  this  approxima- 
tion will  be  discussed  later# 

Collislonal  or  Radiative  Plasma#  Eqs  (2)  or  (3) 
can  be  simplified  if  either  collislonal  or  radiative 
processes  can  be  ignored#  In  a dense  plasma,  collislonal 
processes  should  dominate!  while  radiative  processes  should 
control  the  recombination  in  a rare  plasma#  Zel'dovich 
and  Raiser  (Ref  102 i 408)  derive  an  expression  which  allows 
a plasma  to  be  classified  as  collision  or  radiation  dominant# 
This  derivation  was  accomplished  by  comparing  classically 
derived  rates  for  collisional  and  radiative  recombination# 

Zf 

N#  < 3*1  * 1013  T#3«75  (10) 

v 


13 


for  a hydrogen  plasma,  then  the  radiative  processes  dominate 
compared  to  three-body  collisions*  Thus,  at  low  electron 
temperatures  and  high  electron  densities  collisional 
processes  are  important* 

According  to  an  argument  by  Jones  and  Ali  (Ref  71 i 
3-4),  it  is  impossible  to  produce  population  inversions 
through  radiative  cascade  alone*  The  probability  of  fill- 
ing a level  directly  through  radiation  when  the  electron 
thermal  energy  is  small  compared  to  the  ionisation  energy 
is  proportional  to  (n“1Te"1/2) i while  for  temperatures 
greater  than  the  ionization  energy  it  is  proportional  to 
(n~^T  ”^/2)*  jn  b0th  instances,  the  lower  quantum  levels 
are  more  l.ikely  to  be  filled,  resulting  in  no  population 
inversions* 

The  collisionally-dominated  dense  plasma  is  more 
likely  to  fill  upper  states  then  lower  states  during  its 
recombination*  To  understand  this,  consider  the  principle 
of  detail  balance  which  states  that,  under  equilibrium 
conditions,  the  number  of  particles  leaving  a state  is 
equal  to  the  number  of  particles  entering  that  state 
through  the  inverse  process*  Using  the  notation  of 
Table  1 yields 


Substituting  from  Eq  (9)  for  the  ratio  of  (N,/^2) 
yields 


B(c*,n)  a.  n2  B(n,c*) 


(12) 


The  rate  coefficient  for  collisional  ionization  is 
itself  proportional  to  n2,  producing  this  net  result i 
B(c*,n)  oc  n*  (Ref  681760-763).  The  actual  expressions 
for  the  rate  coefficients  will  be  discussed  in  the  section 
on  rate  coefficients  and  effective  probability  and  will  be 
given  in  Appendix  B,  Rate  Coefficients* 

The  above  results  indicate  that  if  laser  action  is 
desired  from  a recombining  plasma  then  dense  systems  must 
be  used* 
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iv.  solution  Equations 


A computer  program  was  written  to  solve  the  set  of 
equations  which  are  represented  by  Eq  (2).  The  approach 
follows  that  of  Gordiets  et  al  (Refs  36,  37).  The  second 
reference  listed  is  essentially  a translated  condensation 
of  the  Russian  language  paper  which  is  listed  first.  Since 
discrepancies  exist  between  the  two  versions,  comparisons 
of  the  two  works  will  be  facilitated  by  referring  to  them 
as  Paper  I (Ref  36)  and  Paper  II  (Ref  37)  respectively. 

ABBUBrtlQUa 

The  hydrogen  plasma  is  considered  to  be  stationary, 
homogeneous  and  optically  thin  to  radiation,  with  the  ex- 
ception of  the  Lyman  (transitions  to  ground  state)  series 
lines.  The  temperature  of  the  electrons  is  cooled  instan- 
taneously from  an  original  temperature  to  some  final 
temperature.  The  final  temperature  of  the  electrons  is 
assumed  to  be  constant  throughout  the  recombination  process. 
The  use  of  a fixed  temperature  significantly  simplifies 
the  problem.  When  the  temperature  is  fixed,  most  of  the 
rates  need  to  be  calculated  only  once. 

Bq  (2)  will  be  solved  for  the  populations  of  nine 
discrete  levels,  (nj  * 9) • Above  nx,  there  exists  a quasi- 
continuum of  discrete  levels  the  totality  of  which  is 
designated  by  g and  which  will  be  associated  with  the  free 
electron  processes* 
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Imposing  these  limitations*  Eq  (2)  becomes 


dN 

dt 


9 n-1 

-([  £ V(n,m)  ♦ £ R(n,m)]  N 

m=n+l  m=l 


e 


♦ [B(n,c*)  + B(n,g)]  N0  + A*(n)}  Nn 

9 

♦ 2 [R(m,n)N  ♦ A*(m,n)]  N 

m=n+l  8 

n-1 

+ £ V(m,n)N  N 

m»l  e m 

+ [B(c*,n)  ♦ B(g,n)]  Ne3 


♦ [A(c*,n)  + A(g,n)]  N 2 


(13) 


where  n * 1,  2,  • ••*  9* 

The  above  equations  are  supplemented  by  an  equation 
of  conservation  of  total  electrons 


N«  * N_  ♦ N_  ♦ £ N 


g 


n=l 


n 


(14) 


where  Ng  is  the  density  (cm“^)  of  particles  in  the  quasi- 
cnntinuum. 

In  Eq  (13) » A*(n)  is  the  effective  probability  of 
spontaneous  decay  to  any  state  below  n,  whereas  A*(m,n) 
is  the  effective  probability  of  spontaneous  decay  from  the 
state  m to  the  state  n*  The  idea  of  effective  probability 
will  be  discussed  in  the  following  section  on  rate  coeffi- 
cients* 
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Rate  Coefficients  aad  gf£eg.U.Y<?  PrQbflt?i.U$y 

The  coefficients  which  are  necessary  to  solve  Eq  (13) 
are  provided  in  Papers  I and  II  and  are  listed  in  Appendix 
B,  Table  IV. 

The  plasma  which  is  under  consideration  is  assumed 
to  be  optically  thin  to  radiation,  with  the  exception  of 
Lyman  radiation  which  it  may  absorb.  Holstein  (Refs  66, 

6?)  proposes  a modification  of  the  transition  probability 
by  a factor  which  depends  only  upon  the  characteristic 
dimensions  of  the  plasma  and  the  line  shape  of  the  emitted 
radiation  to  explain  the  apparent  trapping  of  radiation 
emitted  from  the  plasma.  The  trapping  of  radiation  within 
the  plasma  is  caused  by  the  existence  of  selective  absorp- 
tion within  the  gas.  A quantum  of  energy  which  finally 
escapes  from  the  plasma  does  so  usually  only  after  having 
undergone  several  absorptions  and  emissions. 

The  effective  probability  is 

A*(m,n)  * F(m,n)  A(m,n)  (15) 

where  F(m,n)  is  the  Holstein  factor.  This  will  be  used 
in  Eq  (13)*  Assume  an  infinitely  long  cylinder  of  plasma, 

R in  radius,  where  only  Lyman  lines  are  affected  and  where 
they  are  Doppler  broadened.  Then, 


F(m,l) 


1 if  k(m,l)R  £ 2 


(16) 


F(m,l) 


if  k(m,l)R  > 2 

k(m,l)R  {nln[k(m,l)R]}1/2  (17) 

where  k(m,l)  is  the  absorption  coefficient. at  the  center 
of  the  spectral  line  (Ref  67 • 1166) • Since  the  absorption 
coefficient , k(m,l),  depends  upon  the  ground  state  popula- 
tion, the  effective  probability  will  also  change  in  time. 

In  Eq  (15),  the  values  of  the  spontaneous  decay 
probability,  A(m,n),  were  obtained  from  the  tables  of 
Weise  si  aJL  (Ref  101»7, 210-211)  when  possible.  The 
appropriate  formula,  Table  IV,  was  used  to  calculate  the 
probability  of  spontaneous  decay  for  the  few  upper  levels 
which  did  not  have  tabulated  values.  Likewise,  the  total 
spontaneous  decay  from  a level  n was  calculated  directly, 
rather  than  by  use  of  the  large  n approximation  (Ref  12 i 
269)  A(n)  * 1.66  x 1010n”4,5  suggested  by  Gordie ts  fit  al. 

Initial  Conditions 

In  all  cases,  the  electrons  were  assumed  to  be  instan- 
taneously cooled  from  a temperature  of  2 eV  to  some  arbi- 
trary lower  electron  temperature.  This  lower  temperature 
was  used  to  determine  the  various  rates*  The  total 
density  chosen  entered  the  problem  as  a parameter.  At  an 
initial  temperature  of  2 eV  and  at  densities  less  than 
lO1^  (cm"^},  the  plasma  is  completely  ionised  (Ref  961 12). 
The  recombination  commenced  at  time  equals  aero  with  all 
the  discrete  levels  vacant.  The  radius  of  the  plasma  was 
assumed  to  be  0*3  cm. 
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Computer  PrflgTHB 

A computer  code  was  written  to  solve  Eq  (13)  and 
the  constraint  equation,  Eq  (14),  for  the  populations 
of  the  nine  discrete  levels  as  a function  of  time*  The 
actual  solution  of  the  differential  equations  was  per- 
formed by  using  a previously  available  program  package* 
This  Runga-Kutta  routine  was  modified  so  that  the  integra- 
tion step  size  was  self-adjusting  during  execution*  See 
Appendix  C,  Computer  Program  Block  Diagram. 


20 


V.  Results  and  Discussion 


Initial  Cgnaitlona 

Computer  solutions  were  generated  for  the  initial 
conditions  listed  in  Table  II* 


TABLE  II 

Initial  Conditions  for  Computer  Solutions 


Ne  (cm"3) 

*.  («V) 

•*# 

lol* 

0.05 

10*5 

0*20 

10* 

0*05 

1016 

0*10 

♦•Calculation  performed  with  and  without 
consideration  of  absorption 

Sain  CalsuIaUgna 

Calculations  of  the  gain  (cm-1)  were  made  for  selected 
lines  at  the  point  and  time  of  maximum  inversion*  The  gain 
was  calculated  from 

0 . 1Q3»10~12  A*(«,n)QwB  - 

A X wr>  (18) 

where  equals  the  wavelength  (cm)  originating  in  state  m 
and  terminating  in  state  m A*(m,n)  is  the  effective  proba- 
bility of  spontaneous  decay  m-»ni  and  equals  the  total 

line  width  of  the  transition  (Ref  3>5)* 

The  line  width  chosen  for  the  calculations  was  the 
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TABLE  III 

Gain  Coefficient  (a)  at  Point  of  Maximum  Inversion 


N#  (cm~3) 

T#  (eV) 

AN' max  (cm*^) 

Transition 

P 

cT 

s, 

h* 

10« 

0.05 

2*2xl014 

CM 

f 

r* 

7*5 

eeio15 

0.05 

14 

2.7xlOl* 

3*2 

9*2 

1015 

0*20 

2.8xl012 

3*2 

•09 

10“ 

0*05 

3*5*1015 

3*2 

25*7 

1016 

0.10 

2.9X1015 

3*2 

21.3 

10“ 

0.10 

3»5xl015 

2-1 

15-8 

tio16 

0.10 

2.0xl015 

2 — 1 

9*0 

**  No  absorption 


t Source i (Ref  74i15).  Gain  is  calculated  for  full  line 
width,  assuming  Stark  broadening  and  no  absorption* 


Stark  width  because  it  is  larger  than  the  Doppler  line 
width  in  the  region  of  interest  and,  consequently,  gives 
a more  pessimistic  value  of  a*  The  Stark  width  was  calcu- 
lated from  (Ref 8 74i7i93«601) 

* 1'*7X10-10  U2-n2)  N.2/3  (19) 

Unlike  the  authors  of  Ref  74  who  take  one  half  of 
Axmn  as  9X1  aPProx^ma‘t^on  t0  U8ed  E<I  Table  III 

reflects  the  full  line  width  in  the  cited  values  of  a» 

On  subsequent  figures,  inversion  density  is  plotted  in 
terms  of 


AN'  « NB  - Njjgg/gyj  (20) 

to  facilitate  the  use  of  Eq  (18).  In  terms  of  the  defini- 
tion provided  for  population  inversion  in  Eq  (1),  Eq  (20) 
may  be  written 

AN'  - gg  AN  (21) 

JlAfifiE  Pain 

Figures  1-4  and  Table  III  indicate  that  sufficient 
gain  is  available  to  make  laser  action  probable*  Addition- 
ally, results  from  the  work  of  Jones  Al  (Ref  74i20)  are 
plotted  for  comparison  in  Fig.  5»  This  latter  work  uses 
different  rate  coefficients  and,  for  the  particular  case 
shown,  ignores  absorption. 

The  reason  for  the  discrepancies  in  the  magnitude  of 


the  inversions  for  the  transitions  shown  and  their  time 
history  is  not  precisely  known*  Probably,  the  use  of 
different  rates  and  the  neglect  of  absorption  accounts  for 
much  of  the  noted  difference*  Unfortunately,  this  author 
did  not  become  aware  of  this  alternate  computer  code  until 
the  latter  stages  of  this  work*  It  would  have  been  informa- 
tive to  use  the  coefficients  of  Jones  a$  a1  within  the 
computer  code  but  time  did  not  permit  this  course  of  action* 
The  rate  coefficients  used  by  Jones  £$  Al  will  be  presented 
within  Appendix  B* 

The  predicted  gain,  although  very  large,  agrees  reason- 
ably well  with  that  of  Jones  jj£  a1*  However,  it  is  orders 
of  magnitude  greater  than  the  predictions  of  the  Russian 
work,  Paper  II*  The  explanation  for  this  discrepancy  is  not 
known,  but  certain  inconsistencies  in  Paper  II  were  de- 
tected* The  Russians  list  gains  for  times  at  which  popula- 
tion inversions  no  longer  existed  according  to  the  graphical 
data  (Ref  37*736-738  and  tabular  data  Ref  36*798-799)* 
Calculations  could  not  reproduce  the  quoted  gain  values  of 
Oordiets  At  Al  even  at  times  at  which  inversions  did  exist* 

The  magnitudes  of  the  gain  estimated  are  very  much 
dependent  upon  the  methods  chosen  to  handle  absorption  and 
cooling* 

Absorption 

Examination  of  Fig*  2 shows  that  absorption  affects 
not  only  the  magnitude  but  also  the  duration  of  the  inver- 
sion between  two  levels*  Suoh  a result  is  certainly 
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- N2gy®2 (c® 


Tine  (nanoseconds) 

Fig*  2 Inversion  Density  of  3— 2 Hydrogen  Transition, 
With  and  Without  Absorption,  vs.  Tine 
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expected  when  the  methods  of  level  depletion  are  considered* 
The  lower  quantum  levels  are  more  likely  to  suffer  radiative 
decay  than  the  upper  levels*  According  to  Bethe,  the  total 
lifetime  of  a level  against  spontaneous  decay  is  propor- 
tional to  n**->  (Ref  12i269)*  It  was  previously  mentioned 
that  the  lower  levels  are  favored  by  radiative  events 
compared  to  the  upper  levels*  Consider  specifically  for  the 
moment  the  3-* 2 transition  (Ha)*  The  lower  level,  2,  is 
decreased  primarily  by  spontaneous  emission*  As  time  pro- 
gresses, the  population  of  the  ground  state,  n = 1,  begins 
to  increase  until  the  effective  probability  of  the  2 -*■  1 
transition  is  affected*  The  population  of  level  2 then 
begins  to  increase  at  a faster  rate  because  the  principal 
avenue  of  depletion  is  reduced*  The  overall  effect,  then, 
is  to  terminate  the  inversion  3 -»•  2 sooner  and  at  a lower 
value  than  would  have  been  reached  had  absorption  been 
ignored* 

However,  the  effect  of  absorption  may,  in  certain 
instances,  increase  the  length  of  duration  of  the  inversion* 
A lengthening  occurs  if  the  upper  level  is  affected  by 
absorption  while  the  lower  level  is  the  ground  state*  If 
inversion  exists  prior  to  the  onset  of  absorption  effects, 
then  the  decrease  in  the  rate  of  upper  level  evacuation  in- 
creases the  rate  of  population  of  level  2*  Simultaneously, 
it  decreases  the  rate  of  growth  of  level  1*  Consequently, 
the  inversion  will  last  longer* 

This  may  partially  explain  some  of  the  discrepancy 
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Fig.  3 Inversion  Density  of  3—2  Hydrogen  Transition 
for  Dlffsrsnt  Electron  Densities  vs.  Tine 
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between  results  of  this  work  and  Ref  74  as  noted  in  the 
comparison  in  Fig.  5. 

It  has  been  assumed  that  stimulated  emission  within 
the  plasma  was  negligible.  This  is  the  approach  taken 
within  the  Russian  literature  in  general.  However,  some 
authors  elsewhere  have  included  this  effect.  Bohn  (Ref  161 
16)  makes  mention  of  these  calculations  but  does  not  publish 
themt  while  Jones  ,§t  al  (Ref  74)  include  stimulated  emission 
explicitly  in  their  computer  code,  although  not  in  the 
instance  cited  in  Fig.  5»  To  the  knowledge  of  this  author, 
there  are  no  presently  available  data  which  allow  for  evalu- 
ation of  the  error  involved  in  ignoring  this  effect. 

Instantaneous 

Although  expedient  for  calculations,  the  assumption 
of  instantaneous  electron  cooling  is  not  experimentally 
justifiable.  Experimentally,  numerous  schemes  have  been 
proposed  to  achieve  the  cooling  of  the  electrons.  These 
include  collision  with  container  walls,  heavy  particle 
collisions,  and  expansions  into  vacuum  and  through  nozzles* 
These  schemes  attempt  to  ensure  that  cooling  is  accomplished 
as  quickly  as  possible.  Since  the  cooling  must  take  place 
in  less  time  than  it  takes  to  produce  the  inversions,  cool- 
ing schemes  may  be  a limiting  factor  in  the  plasma  laser. 

Appendix  D lists  specific  references  within  the  litera- 
ture which  deal  with  the  cooling  problem.  These  indicate 
that  even  with  finite-time-cooling,  inversion  and  subsequent 
gain  are  attainable. 
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An  examination  of  available  literature  indicates  that 
the  Russians  have  a rather  comprehensive  program  for  the 
theoretical  study  and  development  of  the  plasma  laser* 
However,  published  information  from  within  the  U.S.S.R* 
concerning  experimental  development  of  population  inversion 
through  the  use  of  a recombining  plasma  is  surprisingly 
» limited.  Much  of  the  Russian  theoretical  work  has  indicated 

rather  specific  experimental  techniques  (Ref  7»931-93?)» 
but  there  appears  to  be  no  subsequent  experimental  work 
published*  Experimental  confirmation  of  this  principle 
originates  largely  from  outside  Russian  boundaries. 

Hoffman  and  Bohn  (Ref  651878-88O)  have  verified  the 
existence  of  population  inversion  in  expanding  hydrogen 
plasmas  for  the  Paschen  (4-* 3)  and  Brackett  (5—4)  transi- 
tions* The  4 — 3 laser  transition  has  also  been  detected 
in  the  work  of  Bockastein  si  Si  (Ref  15*1260)  working  with 
pulse  discharges  into  a medium  in  which  hydrogen  was  present 
as  an  impurity.  The  line  was  detected  during  the  afterglow 
period. 

The  experimental  demonstration  of  the  feasibility  of 
the  plasma  laser  is  perhaps  best  accomplished  in  the  works 
of  Collins  et  al  (Refs  18il9)«  Working  with  dense  Helium 
plasmas,  they  have  successfully  detected  stimulated  emis- 
sions during  the  afterglow  stage.  These  scientists  maintain 
that  a laser  would  have  been  realized  had  a longer  gain 
path  been  available. 
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VI.  Recommendations 


Presently*  enough  theoretical  and  experimental 
evidence  exists  to  recommend  the  further  study  and  de- 
velopment of  plasma  lasers*  Present  limitations*  such 
as  the  absense  of  rates  for  specific  processes,  will  be 
erased  only  through  experimental  and  developmental  pro- 
grams* 

In  terms  of  the  work  of  this  report,  several  logical 
extensions  suggest  themselves.  Questions  of  plasma  pro- 
duction, electron  cooling  and  material  optimization  which 
were  side-tracked  because  of  time  constraints  certainly 
need  to  be  investigated.  Answers  to  these  questions  will 
help  remove  or  modify  some  of  the  assumptions  made  in 
thiB  study,  Buch  as  the  instantaneous  cooling  of  electrons 
to  a final  arbitrary  temperature. 

In  the  following  pages,  recommendations  are  given 
which  allow  for  extension  of  the  scope  of  this  report. 

The  tasks  outlined  are  neither  necessarily  of  equal  com- 
plexity nor  are  they  of  equal  impact  upon  final  results. 
Rather,  they  represent  what  might  have  been  attempted 
had  more  time  been  available. 

The  computer  code  calculates  the  necessary  rate 
coefficients  of  Eq  (2)  within  the  main  program.  Modi- 
fying the  code  so  that  each  rate  calculation  is  performed 
within  individual  subroutines  would  greatly  facilitate 
the  comparison  of  results  when  different  rate  formulas 
are  employed.  Perhaps  one  of  the  first  attempts  should 
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be  a comparison  of  the  results  obtained  when  using  the 
rate  coefficients  of  Jones  a i (Ref  7*0  versus  the 

rates  of  the  present  study.  By  changing  only  the  rate 
coefficients , the  impact  of  this  change  upon  population 
inversion  may  be  determined* 

Previously,  it  was  mentioned  that  the  assumption  of 
instantaneous  cooling  of  the  electrons  to  some  final 
arbitrary  temperature  could  not  be  defended  on  physical 
grounds.  Possible  schemes  of  electron  cooling  should  be 
investigated.  Expansion  of  the  plasma  or  introduction 
of  heavy  particles  are  possible  methods  of  cooling  the 
electrons*  In  using  a specific  method,  not  only  is  the 
assumption  of  instantaneous  cooling  removed,  but  also  a 
concurrent  time  history  of  the  electron  temperature  is 
provided.  In  terms  of  the  mechanics  of  the  computer 
program,  the  calculation  of  the  various  rates  must  be 
iterated  throughout  the  total  computation  rather  than 
being  calculated  just  once  as  they  are  when  the  final 
electron  temperature  is  constant.  Presently,  the  final 
temperature  of  the  electron  is  arbitrarily  chosen*  The 
use  of  a "physically  realisable"  system  should  remove 
some  of  the  arbitrariness  of  the  electron  temperature  and 
allow  for  laser  operation  predictions  in  terms  of  more 
realistic  parameters. 

The  question  of  using  constituents  other  than 
hydrogen  is  oertainly  not  trivial.  As  noted  on  several 
different  occasions,  even  "well-known"  hydrogen  has  rates 
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that  are  not  understood  exactly*  Attempts  to  scale 
hydrogen  results  to  hydrogen-like  systems  are  frequent- 
ly encountered  within  the  literature  (Refs  10,  11,  16)* 
However,  the  conclusions  drawn  from  these  attempts  are, 
at  best,  likely  to  be  only  as  good  as  the  hydrogen 
results  and  may  be  much  poorer*  By  the  same  token, 
calculations  for  non-hydrogen-like  systems  are  beset 
with  major  difficulties  in  determining  the  appropriate 
rate  coefficients  for  the  various  processes*  Current 
literature  shows  that  the  trend  in  plasma  laser  develop- 
ment is  definitely  towards  these  systems,  but  that  an 
analytic  approach  is  much  more  tedious*  Despite  these 
hindrances,  it  is  probably  in  the  area  of  non-hydrogen- 
like  systems  where  the  plasma  laser  will  finally  be 
constructed* 
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APttgpfliX  A 

Methods  ApprpaipaUon 

Frequently,  Eqs  (2)  and  (3)  must  be  solved  through 
the  use  of  one  approximation  or  more*  Within  this  report, 
specific  techniques  of  approximation  have  been  applied  to 
these  equations*  These  techniques  represent  only  some  of 
those  which  are  cited  in  the  literature  on  the  plasma  laser* 
Other  methods  which  are  not  used  specifically  within  this 
work  but  which  are  often  employed  are  these i one-photon 
approximation,  quasi-stationary  approximation,  and  diffu- 
sion approximation.  These  frequently  used  approaches  will 
be  discussed  briefly  and  references  for  further  study  will 
be  provided*  No  attempt  will  be  made  to  justify  a particu- 
lar method* 

The  works  of  Biberman  sX  al  (Ref  14)  and  Gudzenko  &£ 
al  (Ref  59)  provide  summaries  of  useful  approximations* 

Both  works  list  extensive  bibliographies  for  additional 
information* 

ApregxlBfttifln 

The  form  of  Eqs  (2)  and  (3)  is  simplified  when  the 
one-photon  or  single-quantum  approximation  is  used*  It 
is  predicated  upon  the  fact  that  transitions  are  most  prob- 
able between  adjacent  states*  Consequently,  many  terms 
involving  a given  level  and  another  level  which  is  not 
immediately  adjacent  to  it  may  be  ignored*  Application 
of  the  single-photon  approximation  is  especially  useful 
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when  radiative  processes  can  be  ignored*  Then  this 
approximation  is  used*  only  the  main  diagonal  elements  and 
the  immediately  adjacent  diagonal  elements  in  the  relaxa- 
tion matrix  are  non-zero*  This  allows  the  system  to  be 
solved  in  terms  of  the  populations  of  the  ground  state  and 
free  electrons* 

Referencesi  [13,  14,  59*  68,  98] 

Quasl-stat ionary  ApprqxiTOtiQn 

By  using  the  quasi-stat ionary  approximation,  the 
system  of  rate  equations  can  also  be  algebraically  solved 
in  terms  of  the  ground  state  population  and  the  free  elec- 
tron density.  The  basis  for  this  approach  is  the  fact  that 
under  most  conditions  the  time  for  relaxation  of  the  ground 
state  and  the  free  electrons  is  significantly  longer  than 
the  relaxation  time  of  the  excited  levels  (Ref  10 i 300)* 
Consequently,  all  derivatives  appearing  on  the  left  hand 
side  of  Eqs  (2)  and  (3)  may  be  set  equal  to  zero,  with  the 
exception  of  the  ground  level  rate  which  is  assumed  not  to 
equal  zero*  This  particular  approach  is  employed  very 
frequently  and  is  referred  to  within  the  literature  under 
these  name si  BKM  (Bates,  Kingston,  and  Me Whir ter)  method, 
constant  sink,  stationary  sink,  and  constant  runoff  approxi- 
mation* 

Referencesi  [10,  11,  44,  59,  64,  78] 
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Diffusion  AMfgalfffttlOD 

In  the  diffusion  approximation,  the  upper  levels  are 
characterised  by  quantum  numbers  and  energies  which  are 
assumed  to  be  continuous  functions*  This  approach  is 
based  upon  the  fact  that  at  high  values  of  n the  spacing 
between  levels  gets  smaller  and  smaller*  The  rate  equa- 
tions can  then  be  solved  by  equations  of  the  Fokker-Planck 
type*  Closely  related  to  this  approach  is  the  modified 
diffusion  approximation  (MDA)  which  attempts  to  combine 
the  continuous-like  properties  of  the  highly  excited  states 
with  the  obviously  discrete  nature  of  the  lower  energy 
states  where  the  spacing  between  levels  is  significant*  In 
effect,  it  unites  the  diffusion  and  single  quantum  approxi- 
mations* 

References i £13*  14,  59,  61,  86] 
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Appendix  B 


Rate  ggalf  Is  torts 


The  rate  coefficients  of  Refs  36  and  37  were  used 
to  characterize  the  processes  shown  in  Table  I*  The 
appropriate  formulas  are  listed  in  this  appendix  in  Table 
IV.  The  authors  of  the  above-mentioned  references  attribute 
these  rates  to  semi-empirical  information  and  the  assumption 
of  a Maxwellian  velocity  distribution  of  the  free  electrons 
(Ref  37 » 735) • Unfortunately,  there  is  no  definitive 
information  within  the  available  literature  to  suggest  that 
these  rates  are  more  or  less  correct  than  other  possible 
formulations*  Even  for  the  case  of  hydrogen,  the  appropri- 
ate cross-sect ions  from  which  the  collisional  rates  are 
derived  are  subject  to  question. 

Jones  and  Ali  (Ref  7*0  use  rates  attributable  to 
Drawin  (Ref  27 i 484-486)  for  three-body  recombination  and 
ionization.  Their  three-body  ionization  coefficient  is 
in  the  notation  of  this  report 

B(n,c*)«  1.18xl0"8  n3RyT/*  ^(1,^)  (22) 

where 

*n  " *3/***%  (23) 


and 


*(1*^)  ■ / (l-Xj/uJexpt-u)  In  (1.25uAn)  du  ^ 

*n 
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The  rate  for  the  reverse  process  is  obtained  through 
application  of  the  principle  of  detail  balance*  The  co- 
efficients of  excitation  and  de-excitation  by  electron 
inpact  used  by  Jones  and  Ali  are  attributable  to  Seaton* 
Changing  to  the  notation  of  this  report  yields 

V(n,m)  « 8*2xl0"7  n"^n“^(n“2-m"2)"1*'  T^expMJ^-XJ  ) ^5) 

while  R(m,n)  comes  from  detail  balance  considerations* 

Numerical  evaluation  of  the  coefficients  of  Gordiets 
et  al  compared  to  those  employed  by  Jones  indicate  that 
the  former  values  are - generally  smaller  by  a factor  of 
102  in  the  region  of  interest.  This  discrepancy  is 
apparently  significant!  but  in  terms  of  the  usage  in  this 
report*  it  indicates  that  the  very  existence  of  an  inversion 
is  not  very  dependent  upon  the  rates  used*  However*  the 
impact  on  gain  and  duration  of  inversion  may  be  significant 
for  different  rates  (Ref  73).  Obviously,  experimental 
verification  of  the  appropriate  rates  in  dense  plasmas  is 
needed* 

The  rates  of  Gordiets  at  al  were  used  for  two  reasons! 
First*  there  appeared  to  be  no  priori  evidence  contradicting 
their  use  for  estimate  purposes!  and,  second,  possible  con- 
firmation of  a numerical  solution  might  be  obtained  by 
standardisation  of  rates* 


^9 


TABUS  IV 

8 For  Rate  Coefficients  (Gordiets  al) 


Sources i (Ref  361 792-79*0  and  (Ref  37»73^ 
Notes  Last  two  rates  are  incorrect  in  Ref 


Pig.  6 Computer  Program  Block  Diagram 
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Bibliography  STQBB-RflXfirgnfig 

The  following  list  of  plasma  laser  interest  areas 
is  cross-referenced  to  the  bibliography  by  number*  The 
topic  under  which  a given  reference  is  listed  is  the 
emphasis  area  of  the  paper*  In  some  cases • the  same 
reference  may  be  listed  under  more  than  one  subject*  Many 
of  the  references  listed  under  "Rate  Coefficients"  do  not 
deal  directly  with  plasma  laser  systems,  but  they  do  pro- 
vide needed  background* 

1.  6,  23.  34.  38.  57. 
75.  79. 

34,  59.  61,  86,  98* 


Alkali  Plasmas  for  Lasers i 
Approximations i 


Chemical  Plasma  Systems i 


33  < 


Cooling  of  Plasmasi 


35.  50.  54.  75.  84,  94. 


Experimental  Worki 


7.  15.  17.  18,  19.  20, 
21,  52.  65,  80,  87,  90* 


Gain  Projections! 


9.  16,  22,  36,  37,  38, 
40.  41,  42.  51.  72,  74. 
91.  95* 


Laser  Devices! 


18.  19.  49.  55.  83,  84. 
85.  91. 


Non-Atomic  Recombination  Principle!  31,  33,  37.  44,  45,  46, 

47.  48,  55.  56,  60.  89* 
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Plasma  Generations 


7.  55.  95* 
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